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Abstract 

Background: Our interest in Candida albicans nnitochondria began with the identification of GOAL We 
demonstrated its role in cell energy production, cross-talk among mitochondria and peroxisomes, non-glucose 
energy metabolism, maintenance of stationary phase growth, and prevention of premature apoptosis. Its absence 
results in avirulence. However, what regulated transcription of GOAl was unknown. 

Results: To identify transcriptional regulators (TRs) of GOAl, we screened a C albicans TF knockout library (TRKO) 
and identified Rbfl p, Hfll p, and Dpb4p as positive TRs of GOAl. The phenotypes of each mutant (reduced 
respiration, inability to grow on glycerol, reduced ETC CI and CIV activities) are reasonable evidence for their 
required roles especially in mitochondrial functions. While the integration of mitochondria with cell metabolic 
activities is presumed to occur, there is minimal information on this subject at the genome level. Therefore, 
microarray analysis was used to provide this information for each TR mutant. Transcriptional profiles of Rbflp and 
Hfllp are more similar than that of Dpn4p. Our data demonstrate common and also gene-specific regulatory 
functions for each TR. We establish their roles in carbon metabolism, stress adaptation, cell wall synthesis, 
transporter efflux, peroxisomal metabolism, phospholipid synthesis, rRNA processing, and nuclear/mtDNA 
replication. 

Conclusions: The TRs regulate a number of common genes but each also regulates specific gene transcription. 
These data for the first time create a genome roadmap that can be used to integrate mitochondria with other cell 
processes. Of interest, the TRs are fungal-specific, warranting consideration as antifungal drug targets. 

Keywords: Transcription factor. Non-glucose carbon metabolism. Mitochondria, Lipid oxidation. Metabolic 
regulation, Candida albicans 



Background past few decades primarily due to cancer chemotherapy, 

Fungal invasive infections of humans are now referred organ/bone marrow transplantation, surgical interven- 

to as "hidden killers" [1]. More than 90% of these infec- tion, and the AIDS pandemic [3,4]. Treatment of these 

tions are caused by species of Candida, Cryptococcus, infections costs ~ $2.0-2,6 billion per year [5-8]. Global 

Aspergillus, and Pneumocystis [1]. Blood-borne, nosoco- cryptococcal meningitis (usually caused by C. neofor- 

mial candidiasis is ranked 4^^ in frequency in the USA mans) in HIV/ AIDS patients and others with immuno- 

with a crude and attributable mortality of 49 and 27% suppression therapies is estimated at 1 million cases per 

(USA), similar to other developed countries [1,2]. The year; 620,000 deaths alone are in Sub-Saharan Africa [1]. 

incidence of candidiasis has increased sharply over the Cryptococcus gattii is an emerging pathogen of appar- 
ently healthy people, signif)^ing its potential as an even 

— \ — \ more dangerous invasive fungus. Death from the top 10 
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that of tuberculosis (1.4 million) or malaria (971,000) 
[1]. The enormity of fungal infections is magnified by 
the non-invasive (superficial) infections such as nail and 
scalp infections (1.9 billion), vaginal infections of women 
during child-bearing years (frequency of 50-75%), and 
oral and esophageal candidiasis in HIV/ AIDS patients 
(12 million). In part, the rising costs are associated with 
inappropriate therapy, defined as delayed intervention, 
inadequate dosage, or administration of an antifungal to 
which an isolate was considered drug resistant [6]. 

C. albicans remains as the most common cause of 
candidiasis among all Candidia species. Virulence of this 
organism is commonly attributed to factors that initiate 
colonization of host cells (the ALS gene family and 
others), cause invasion (secreted lipases and proteases), 
regulate morphogenesis (the yeast hyphal transition), 
and biofilm formation [1]. In vivo virulence of these fac- 
tors has been established in animal models fulfilling the 
paradigm of "Molecular Kochs postulates". Aside from 
the construction of single mutants to verif)^ a role in 
pathogenesis, another useful approach to understanding 
virulence is to characterize global gene differences be- 
tween a pathogen (C. albicans) and a non-pathogen 
(Saccharomyces cerevisiae, model yeast) or between two 
pathogens, one with a much lower incidence of causing 
candidiasis (C. dublinensis) [9]. Both types of data sug- 
gest interpretations of the gene repertoire needed by a 
pathogen. One of the major differences between C. 
albicans and model yeast is a rewiring of transcriptional 
regulation [10]. For C. albicans^ enzymes of alternative 
carbon metabolism (non-glucose substrates) are stabi- 
lized even in the presence of glucose, compared to 
model yeast of which these same enzymes are regulated 
by glucose-repressible events [11]. Speculation is that C. 
albicans maintains a backup source for energy and 
carbon conservation to respire when confronted with 
low levels of host glucose. Model yeast when grown 
aerobically uses glucose via glycolysis and is referred 
to as Crabtree-positive. Oppositely, C. albicans respires 
oxidatively in the presence of glucose and is Crabtree- 
negative [12]. These observations are not surprising, 
given the differences in their environmental niches. In 
the case of C. albicans, low blood levels of glucose cause 
the utilization of alternative carbon sources as men- 
tioned above and described in other labs [13-16]. Some 
peroxisomal activities in C. albicans are critical to the 
pathogenesis of candidiasis, since these organelles house 
alternative carbon metabolic pathways (such as the 
glyoxylate cycle) that are critical to survival of the organ- 
isms in macrophages [15]. 

Our interest in mitochondria of C. albicans began with 
the identification of GOAl [16]. Functional annotation 
was developed based upon phenotypic assays of a goalt^ 
null mutant. Goalp translocates to mitochondria during 



stress and in the presence of non-glucose substrates 
such as glycerol. The protein regulates complex I (CI) of 
the electron transport chain and also interacts with 
peroxisomes, organelles that house alternative carbon 
metabolic pathways. The loss of GOAl causes a major 
reduction in mitochondrial membrane potential and a 
concomitant reduction in the formation of ATP. We 
have shown that a dysfunctional CI causes an increase in 
reactive oxidant species (ROS), triggering apoptosis and 
an associated shortened chronological aging in vitro 
[16-18]. We demonstrated that cross-talking between 
mitochondria and peroxisomes in the presence of either 
glucose or non-glucose substrates requires Goalp [17,19]. 
Importantly, there are several defects in the mutant in re- 
gard to virulence and host cell interactions [20]. Compared 
to parental and gene-reconstituted strains, goal^ is aviru- 
lent in a murine model of blood-borne candidiasis, killed 
readily by human neutrophils ex vivo, and hypersusceptible 
to triazoles. Microarray analyses of goal^ indicate a major 
down regulation of genes associated with peroxisomal func- 
tions and carbon metabolism [17,19]. But, what regulates 
GOAP, Herein, we report on the identification of three 
transcription regulator (TR) Rbflp, Hfllp, and Dpb4p, each 
of which positively regulates GOAl, Importantly, we assign 
functions to each of these transcription factors in the regu- 
lation of cellular processes. Two of these TRs have not been 
functionally annotated, and the third (Rbflp) is a repressor 
of filamentation [21,22]. Using microarray analysis, we 
show that there is both TR-specific gene regulation as well 
as regulation of a set of common genes. 

Results 

Identification of transcription regulators of G0A1 

To initially identify transcription factors that regulate 
GOAl, we screened a transcription regulator knockout 
(TRKO) library of C. albicans [21]. The library was 
maintained in 96-well microtiter plates at -80°C. The 
initial screen of 163 TR mutants was done in 36-well 
plates containing YP-dextrose (YPD) or YP-glycerol 
(YPG), since glycerol is only a substrate for mitochon- 
drial oxidation. Cultures of each mutant were grown 
overnight and inoculated in each medium. We identi- 
fied 6 mutants that were either unable to grow or 
grew poorly only on YP-glycerol (YPG) (Figure lA and 
Table 1). Two independently generated mutants of 
each gene were used in assays for growth on both YPD 
and YPG (Figure lA). Each independent gene mutant 
is identified as X and Y. 

The mutants that did not grow on YPG were assayed 
for transcription of GOAl using real-time PCR (Figure IB). 
Of the six mutants, we observed significant down regula- 
tion of GOAl in mutants lacking HFLl, RBFl, or DPB4, 
Transcription levels for both GOAl and NDH51 (Ndh51p 
is a complex I subunit protein) were decreased 4-6 fold in 
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Figure 1 Glycerol utilization and transcriptional regulation of G0A1 and NDH51 by the TRKOs. A. Drop plate assays of 6 TRKO mutants in 
YP-dextrose (YPD) or YP-glycerol (YPG). Each mutant was spot inoculated with 3 \i\ of stock cells and grown at 30°C for 24 h. X- and Y-sets refer 
to two independently constructed mutants lacking each gene. The reduction or lack of growth on YPG indicates a mitochondrial respiratory 
dysfunction. B. Each of the mutants shown in (A) was measured for transcription of GO A 1 by Real-time PGR. Data are indicated as a down 
regulation of either GOAl or NDH51. The latter encodes a Gomplex I subunit of the mitochondrial electron transport chain. Dpb4p, Hfllp and 
Rbflp are positive regulators of both GOAl and NDH51. 



RBFl, HFLl and DPB4 TRKO strains. Orfl9.2088 is 
named DPB4 which corresponds to the S. cerevisiae ortho- 
log. The MACl and HCMl mutants had much smaller 
changes than the RBFl, HFLl, or DPB4 TRKO strains vs. 
wild type (WT) cells (Figure IB). The knockout strain for 
TUPl was excluded from further studies because of its 



poor viability. HCMl and MACl are conserved transcrip- 
tional regulators in C. albicans and S, cerevisiae [21]. 
Maclp and Tuplp are thought to regulate copper and iron 
uptake [23]. Our data indicate that we have identified three 
TRs that regulate GOAl and NDH51 both of which are 
required for mitochondrial activity in C. albicans. Of 



Table 1 TRKOs those are unable to grow on YP-Glycerol and their mutant phenotypes 



ORF 


Gene 


Description 


Mutant phenotype 


orfl 9.4853 


HCMl 


Has a forkhead domain; involved in hyphal morphogenesis; similar to the S.cerevisioe 
ortholog that regulates transcription during S-phase of the mitotic cell cycle 


Vegetative growth decreased in YPG 


orf 19.7068 


MACl 


Regulon of CoCTR], a high affinity copper transporter gene; role in ion homeostasis 


Vegetative growth decreased in YPG 


orfl 9.2088 


DPB4 


Putative DNA polymerase epsilon subunit D; Ortholog of S.cerevisioe involved in 
DNA-dependent DNA replication 


Viable, slow-growing and displays 
abnormal invasive growth 


orfl 9.3063 


HFLl 


HAP5-like; phosphorylated protein 


Vegetative growth decreased in YPG 


orfl 9.5558 


RBFl 


Glutamine-rich activation domain; binds RPG-box DNA sequences; antigenic during 
human oral infection 


Vegetative growth decreased in YPG 


orfl 9.6 109 


TUPl 


Co-repressor; represses filamentous growth; regulates switching 


Abnormal colony appearance on 
YPD and YPG 
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Figure 2 The yeast to hyphal transition, generation time, and transporter activity are altered in each of the TRKO mutants. A. Growth 
on Spider agar, YPD pH 9.5 and YPD pH 5.0. On Spider agar, all mutants had reduced colony growth. On YPD (below), colony diameter was again 
reduced and each mutant had a wrinkled colony appearance with a lack of filamentation at pH 9.5. B. Microscopic growth of each mutant is 
shown. In YPD (30°C), each of the TRKO mutants was constitutively filamentous compared to the yeast growth of WT cells (SN250) as reported 
previously [21,22]. In 10% serum, all mutants and WT strains are filamentous. C. Doubling times are shown for WT (SN250) as well as each TRKO 
strain. Strains dpb4I\ and rbflA were significantly higher than WT cells. Strain hfllA was not significantly different than WT cells. Doubling time was 
determined by cell mass measurements. 
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importance, their role as regulators of cell metabolism has 
not been described The remaining sections reflect our 
studies of the 3 TR mutants whose transcription of GOAl 
was down regulated. 

Morphology and growth of hfllA, rbfiA, and dpb4A 

A comparison of morphology as weU as doubling times 
was carried out with the RBFl, HFLl and DPB4 TRKO 
strains (Figure 2,A-C and Table 1). Abnormal colony 
phenotypes (Figure 2A) and filamentous microscopic 
growth (Figure 2B) were observed in most growth condi- 
tions compared to the parental strain. Each of the three 
TRKO strains exhibited smaUer colonies and a near ab- 
sence of filamentous rings at the edges of colonies com- 
pared to WT colonies on Spider medium (Figure 2A). 
The surface of colonies from mutant strains was extensively 
wrinkled by day 7 on both Spider and YPD (pH 9.5) agar 
media. Invasive growth on YPD-2% glucose, pH 9.5, also 
was diminished in the three TRKO strains. All mutants 
were constitutively filamentous in YPD at 30°C. However, 
dpb4A was less filamentous compared to the other two 
mutants (50% of cells), and those were mostly pseudohy- 
phae. RbflA and hfllA were similar to WT cells in 10% 
serum at 37°C, as reported previously for rbflA [21]. The 
growth phenotypes of each mutant are summarized in 
Table 1 and also described at the Candida CGD database 
(www.candidagenome.org). Doubling times varied accord- 
ing to the specific mutant from 3.32 to 5.32 hr compared to 
WT cells (2.84 hours) (Figure 2C). 

Antifungal susceptibilities of the TR mutants 

Since we showed previously that deletions of GOAl and 
NDHSl resulted in hypersensitivity to triazole antifun- 
gals [19], we compared the susceptibflities of each TR 
mutant to the antifungals fiuconazole, amphotericin B 
(AmB), and caspofungin (Table 2). Since trafling has 
been reported, we measured both MIC50 and MICioo for 
fluconazole. Simflar to our published data on goal A 
and ndhSlA, rbflA and hfllA are hypersusceptible to 
1-2.0 (ig/ml fluconazole (MICioo)> while dpb4A was 
similarly susceptible as parental cells. The susceptibilities 

Table 2 Antifungal drug susceptibilities to cell membrane 
and cell wall inhibitors^ 



Gene 
KO 


Fluconazole(|jg/ml) 
MICioo IVIIC50 


AmB (ng/ml) 
MICioo 


Caspofungin (ng/ml) 
MICioo 


WT 


256 


0.5 


0.25 


0.25 


GOAl 


1.0 


1.0 


1.0 


0.0625 


DPB4 


256 


1.0 


1.0 


0.5 


HFLl 


2.0 


0.25 


0.5 


0.125 


RBFl 


2.0 


0.25 


0.25 


0.0625 



of aU mutants to AmpB and caspofungin were not statis- 
tically different compared to parental ceUs. 

The TF mutants are hypersusceptible to inhibitors of cell 
wall formation 

Surprisingly, all mutants were hypersusceptible to calco- 
fluor white (CFW) and caspofungin on drop plate assays, 
but rbflA^LYid hfllA were more so than dpb4A (Figure 3). 
As for susceptibility to Congo red (a p-l,3-glucan inhibitor), 
only the rbflA and to a lesser extent hfllA were more sus- 
ceptible than WT cells. Thus, rbflA was affected most by 
cell wall formation inhibitors. The differences between the 
caspofiangin MICs (described above) and cell wall inhibitor 
drop plate assays suggest that the regulation of cell 
wall integrity among the three TRs is different. How- 
ever, the disparities of both assays could also be an 
explanation as MIC determinations were done using 
RPMI medium at 37°C whUe drop plate assays were 
done in YPD agar at 30°C. 

Mitochondrial defects in rbflA, hfllA and dpb4A 

The inability of each mutant to assimilate a non- 
fermentable carbon source (glycerol) indicates defects in 
mitochondrial respiration. Further, we were interested in 
comparing the functions of each of the TRKO strains to 
GOAl in energy production and carbon metabolism. To 
determine the mitochondrial status of the TRKOs, we 
first measured oxygen consumption among mutants and 
parental cefls. The oxygen consumption rates were de- 
creased by 2.2-fold for dpb4A, and about 5-fold in hfllA 
and rbflA compared to WT cells (Figure 4A). For these 



*Data reflect averages of 3 experiments. 
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Figure 3 The TRKO mutants (indicated on the left) are 
susceptible to the cell wall inhibitors calcofluor white, Congo 
red and caspofungin. Cell wall inhibitors are YPD only (control), 
YPD + caspofungin (CsP) first row; YPD + calcofluor white (CFW) or + 
Congo Red (CR) in the second row. Compared to growth on YPD 
only, each mutant was hypersensitive to CFW and caspofungin, 
while only r6f/A was hypersensitive to CR. The hf//A displayed a 
slight sensitivity to CR. 
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Figure 4 The TRKO mutants are deficient in respiratory functions. A. respiration; B. ETC CI activity; C. ETC CIV activity, and D. ROS production. 
dpMA, hfllA, and rbflA each respired less (2.2-5-fold reduction) than WT cells. CI activity decreased in each mutant proportionally to their 
decrease in respiration. CIV activity was 5.5 fold lower in the TRKO mutants. In D, ROS production compared to WT cells was highest in the rbflA, 
although also significantly increased in hfllL The dpb4/\ produced ROS equal to WT cells. E. ETC CI and CIV inhibitors: growth of mutants and WT 
cells on YPD containing either rotenone (CI inhibitor) or KCN (CIV inhibitor) is shown. rbflA and hfllA are hypersensitive to both inhibitors while 
dpb4/\ was less so. 



experiments, total oxygen consumption was determined 
from equal masses of cells (per mg dry cell mass, DCM). 

The ETC CI and CIV activities (Figure 4B, C), reactive 
oxidant levels (ROS) (Figure 4D, and susceptibilities to 
ETC CI and CIV inhibitors (Figure 4E) were also 
evaluated in rbflA, hfllA and dpb4A compared to WT 
cells. ETC CI and CIV enzyme activities for the rbflA 
mutant were significantly reduced by 4-fold and 14-fold, 
respectively. Corresponding to the decrease in CI en- 
zyme activity was an increase in sensitivity to rotenone, 
a CI inhibitor and KCN (CIV inhibitor) in rbflA, For 
hfllA, CI activity was less affected than rbflA, but CIV 
activity was reduced similarly to rbflA. CI enzyme activ- 
ity in dpb4A was similar to that of hfllA. Sensitivity of 
the dpb4A to rotenone was less than that of the other 2 



mutants but the same as hfllA in regard to KCN sensitivity. 
These data indicate that each of the TR mutants have 
altered CI and especially CIV enzyme activity although 
correlates with complex inhibitors are not absolute. 

One of the striking features of mitochondria with 
dysfunctional CI and CIV activities of the ETC is an in- 
crease in mitochondrial ROS [17,18]. In this regard, ROS 
levels were nearly 20-fold higher in rbflA and ~5-fold 
higher in hfllA; however, ROS production in dpb4A was 
similar to that of parental cells (Figure 4D), indicating 
that the ROS scavenging system was less functional in 
hfllA and rbflA but not affected in dpb4A, Microarray 
data indicated that genes associated with ROS detoxifi- 
cation such as SODS, GPXl, GPX2, in each mutant were 
increased slightly, but a down regulation in S0D6 and 



Khamooshi et al. BMC Genomics 2014, 15:56 
http://www.bionnedcentral.conn/1471 -21 64/1 5/56 



Page 7 of 18 



GRXl occurred in both hfllA and rbflA (Additional file 1: 
Table SI, Additional file 2: Table S2 and Additional file 3: 
Table S3). The decrease in S0D6 and GRXl transcrip- 
tion may partially explain the high ROS levels in hfllA 
and rbflA, 

Global transcriptional profiling in rbfiA, hfllA, and dpb4l^ 

Based upon our published data on transcriptional profiling 
of the goal^ [19] and the functions of the RBFl, HFLl, and 
DPB4 as positive regulators of GOAl, we expected com- 
mon gene pools as well as TR-specific gene changes. To ob- 
tain data to support this premise, we compared array data 
from each TR mutant to goal A versus their own parental 
strains. A 2-fold increase/decrease in transcription was used 
to determine if significant changes occurred. 



General observations of transcriptional changes for 
each TR mutant 

The total number of genes whose transcription changed 
significantly compared to SN250 was 862 (rbflA), 692 
{hfllA) and 505 {dpb4A) (Figure 5). The genes with up/ 
down changes in expression vs, the parental strain were 
grouped for each TR mutant based on their functional 
classification. The assignment of functional categories is 
based on the information provided by the C. albicans 
CGD and S, cerevisiae databases. The% of uncharacter- 
ized genes in hfllA and rbflA was 27-28%, while that of 
dpb4A was 17%. We found that 20% of the genes from 
rbflA and hfllA were associated with mitochondria, car- 
bon, lipid and amino acid metabolic process. For dpb4A, 
the genes of metabolic processes accounted for only 14% 
(Figure 5). 





Figure 5 Total gene changes in the TRKO mutants of C. albicans rbfll^ (A), M7A (B) dpb4l^ (C). Genes with more than a 2-fold or less than 
2-fold change in expression levels are included in the diagram. Functional classification of each category is based on the C. albicans genome 
database (http://www.candidagenome.org/) and curated manually with 5. cerevisiae genome database (http://www.yeastgenome.org/) 
when needed. 

V J 
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To determine differences among the three TRs in 
metabolic genes, we compared each TR transcriptome 
with the goal A transcription profile. We found that a 
significant overlap of common genes was observed be- 
tween rbflA and goal A (334 in total) and between hfllA 
and goal A (302 in total). By contrast, the number of 
genes common to both dpb4A and goalA was less (247 
in total). Gene data are represented in separate sections 
in Figure 6A-C for each TR mutant in comparison to 
goal A data published previously [19], and Figure 6D 
compares rbflA and hfllA, Common down regulated 
genes encoding putative metabolism functions cluster 
in the lower left quadrants for the three TR mutants 
(Figure 6A-C), respectively. Common cell wall, morphology 



switching, and stress upregulated genes cluster in the upper 
left quadrants. The group of genes in the lower right quad- 
rant in dpb4A (Figure 6C) is related to ribosomal /mtRNA 
processing and DNA/mtDNA replication or maintenance, 
which is down regulated in dpb4A, The upper right quad- 
rants for each TR mutant represent a number of transcrip- 
tionally altered but non-functionally clustered genes. Green 
triangles indicate TR-specific genes and red squares in- 
dicate G0A7-specific genes (Figure 6A-C). The down 
regulated cluster of genes in dpb4A only (Y axis of 
Figure 6C) includes 5 mtDNA genes that encode the 
ETC CI subunit. The similarity of gene transcriptome 
changes between rbflA and hfllA reaches a maximum 
in the co-linearized rate (R^ = 0.76) that is derived 
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Figure 6 A comparison of tiie transcriptomes reveal differences among tlie 3 TRKO mutants. Data are derived from RBFl (A), HFLl (B), 
DPB4 (C) TRKO strains vs. the GOAl mutant, and rbfl vs. hfll mutants (D). Genes common to botli tine gool/\ and eitlier rbfl/\, hfn/\, or dpb4/\ are 
depicted as blue diamonds. Tine goa/Z\-specific genes are indicated as red squares, and TRKO-specific genes as green triangles. Upregulated genes 
for each TR mutant are show above the horizontal axis, down regulated genes below the horizontal axis; while the right side and left side of 
vertical axis represent upregulated and down regulated genes for goolI\. Each quadrant reflects a cluster of similar genes between goal A and TR 
mutants. Upper left: cell wall, stress adaptation, and morphological switching; lower left: metabolism; lower right: different rewiring responses 
for common genes, in Figure 6C reflecting DNA replication or maintenance and rRNA processing between dpb4/\ and goal A, Upper right: 
upregulated genes of different functional annotation. 
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from common gene pools in a total of 326 genes 
(Figure 6D). 

There are several points of interest from these data. 
First, down regulated metabolic genes (left lower quad- 
rants) are much greater in number from rbflA (Figure 6 A) 
and hfllA (Figure 6B) than in the dpb^A (Figure 6C). Sec- 
ond, dpb^A has a much greater number of down regulated 
genes that are associated with mitochondria DNA and pro- 
tein synthesis machinery, such as mtDNA encoding genes 
and mitochondrial rRNA ribosomal assembly (lower right 
quadrant) compared to the other TRs (Figure 6C, Table 3). 
As shown in Figure 5, the total% of gene changes in rRNA 
processing and ribosomal assembly was 31% in dpb4A 
compared to 3% in rbflA and hfllA, Third, the cluster size 
of upregulated genes associated with cell wall integrity, 
stress responses, morphological switching, and adhesion 
are reasonably similar for all TRs but opposite with goal A 
(upper left quadrant). These data demonstrate that the 



number of genes associated with metabolic and other cellu- 
lar bioprocesses varies among the TR mutants, suggesting 
differences in their regulatory roles. 

Functional clustering of gene changes In rbfiA, hfllA, 
and dpb4L 

Mitochondrial respiration 

As shown in Figure 4A, oxygen consumption was re- 
duced by 5 -fold of WT levels in the hfllA and rbflA and 
by 2.2-fold in dpb4A, This observation is consistent with 
their increased doubling times in vitro in Figure 2C, 
while also indicating some degree of dysfunctional mito- 
chondria. To better understand the reasons for the 
reduction in respiration, transcriptional profiling was 
performed. Statistically significant changes in mitochon- 
dria genes are summarized among the three mutants 
(Table 3). 



Table 3 The transcription profiles of mitochondrial functions of TRKOs 



Biological processes 




rbflA 


hfllA 




dpb4A 






down 


up 


down 


up 


down 


up 


ETC subunits & assembly, alternative oxidase 


Offl9.1691 


NAD5, A0X2 


NDEl 


NAD6, 
MRFl 


NAD2, NAD3 


MRFl 




ATP 18^ CYTl 


orf 19.3353 




C0X2 


NAD4L, NAD4 






OCRS, 0CR9 


orf 19.6374 




orfl 9.230 


NAD5, ATP 19 






RIP I 
Off 19.2249 


MRFl 




orfl 9.3353 


orfl 9.4273 




Mitochondrial genome maintenance 


SUV3 


PRP13,0GG1 
RFC4 


SUV3,CPA1 
orfl9.3062,RAD32 




SUV3, RPPl 




Mitochondrial tRNA processing & ribosome 
assembly 


orfl9.1201 
orfl9.1356 
orfl9.5161 


MRPL6 


orfl 9.1 356 




orfl 9.2309, 
orfl 9.2384 

orfl 9.3303, 
orfl 9.41 60 

orfl 9.41 4, TRMl 




Mitochondrial carbon metabolism 


QCV2, ALD6 


orfl 9.850 


UCFhCTNl 

orfl 9.2 184, ALD6 

ILV6,DLD1 DLD2, 
PDK2 


orfl 9.3043 
GDB1,MAL2 


IDH2 




Mitochondrial transporters 


Off 19.4583 
orfl9.5129 
orf 19.6555 


PAM17, ERVl 


ATM lorf 19.28 

orfl 9.6532 

orfl 9.6555 
orfl 9.4966 

MDL2,FRP5 

orfl 9.3455 




RIM2, PET9 
orfl 9.5628 
orfl 9.28 




Others 


ISA2 
Off 19.2749 
orf 19.6550 

orf 19.5270 


orfl 9.2401 

orfl 9.2544 

orfl 9.2825, 
orfl 9.4523 

orfl 9.461 2, 
orfl 9.6283 


orfl 9.21 7 5 
orfl 9.6736 
GTT12 


PST2 
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Although hfllA and rbflA displayed the same reduc- 
tion in oxygen consumption, mechanistically the de- 
crease in respiration was not identical For example, 
the genes that are required for assembly of ETC CIV 
{orfl9.1691) and CV {ATP18, orfl9.2249) were down 
regulated only in rbflA (Table 3). Other respiratory-re- 
lated genes such as the ubiquinol cytochrome c reduc- 
tase {QCR8 and QCR9) were down-regulated in rbflA as 
well Candida species are known to have at least 3 re- 
spiratory pathways in addition to the classical ETC chain 
[18,24]. The other pathways (AOX and PAR) are utilized 
when the classical pathway is not functioning well. Thus, 
in rbflls., an impairment of mitochondrial ETC chain 
results in subtle increases of NADS, orf 19.3353 and a 
significant increase of A0X2 (11. 5 -fold). However, these 
apparent compensatory changes in the ETC and alterna- 
tive respiratory route are not sufficient to overcome the 
mitochondrial dysfunctions in rbflA, This is especially true 
for A0X2, which is a hallmark gene of a mitochondrial 
stress response rather than its role in respiration [25]. 

The transcription profile of genes related to mito- 
chondrial respiration in hfllA was different. Instead of 
a direct impact on ETC complex subunits, we found 
that the down-regulated genes in this mutant were 
more related to mitochondrial carbon metabolism such 
as pyruvate {ALD6, PDK2\ lactate {DLDl, DLD2) and 
amino acids (orfl9.2184, ILV6) (Table 3); also four mito- 
chondrial mtDNA maintenance genes (orfl9.3062, SUV3, 
CPAl, RAD32) were transcriptionally reduced in hfllA, At 
same time, eight mitochondrial transporters of amino acid/ 
protein, ammonium, FAD, and magnesium are down regu- 
lated in hfllA, Similar to the rbflA mutant, a possible com- 
pensatory response of mitochondrial respiration in hfllA 
was reflected by up-regulation of CI {NAD6), CIV {C0X4) 
and others (Table 3). For this same mutant, the only down 
regulated gene with a respiratory function was NDEl 
(NADH dehydrogenase), which functions as an ETC CI 
subunit in S, cerevisiae. 

In contrast, the absence of DPB4 suppressed expression 
oiNAD2, NADS, NAD4L, NAD4, and NADS more than 2- 
fold; these five genes comprise the mitochondrial respir- 
ation chain CI which are encoded by mitochondrial DNA. 
Also the mitochondrial nucleotide transporter genes 
(orfl9.28, and RIM2) are down regulated in dpb4A (Table 3). 
Meanwhile, the house-keeping genes for mtDNA mainten- 
ance and mitochondrial RNA processing (a total of 8 genes) 
were transcriptionally reduced (Table 3). Obviously, Dpb4p 
regulates respiration in a different manner than the other 
two TRs. 

Fatty add cataboHsm is correlated transcriptionally with 
decreases in phospholipid (PL) biosynthetic encoding genes 

Simflar to mammalian ceUs, in C. albicans lipids provide 
a source for energy generation via catabolism as well as 



phospholipid biosynthesis via anabolic pathways [26]. 
Regulation of both catabolic and anabolic pathways is 
critical to cell growth [27]. After comparing the tran- 
scriptome of lipid metabolism with goal A, differences 
are seen among the three TR mutants of C. albicans. 
The absence of DPB4 resulted in an upregulation of |3- 
oxidation (lipid catabolism) and genes of the peroxi- 
somal glyoxylate cycle (Table 4). But its PL biosynthesis 
may be compromised since IN04 (PL biogenesis activa- 
tor) was down regulated by 100-fold vs. WT ceUs. The 
other TRKO strains {rbflA and hfllA) resembled goal A, 
and each other, with significant down regulation in lipid 
oxidation, lipase, the glyoxylate cycle, and peroxisomal 
importing systems such as the peroxins. In addition, 
genes for PL biosynthesis including sphingolipid (SL) 
biosynthesis were down regulated whfle genes for PL 
catabolic processes were up regulated. In contrast to the 
DPB4 mutant that may regulate PL biosynthetic process, 
decreased gene expression for lipid catabolism and PL 
biosynthesis in the other two mutants indicate that 
RBFl and HFLl positively regulate both lipid catabolism 
and PL biosynthesis. 

Alternative carbon source metabolism is also regulated 
by each TR 

The biological implications for the assimflation of non- 
glucose carbon sources even when glucose is not limit- 
ing for C. albicans has been described [12,28-30]. We 
observed that numerous genes, required for non-glucose 
utilization in both rbflA (26 of a total of 31 genes) and 
hfllA (23 of 32), were down regulated along with mito- 
chondrial defects. Notably, the GAL gene cluster was 
significantly reduced by 4.6-6.4 fold in hfllA {GALl, 7, 
10, 102) and 2.9-3.0-fold in rbflA {GALl 10) (Table 4). 
On the other hand, most of the genes for alternative car- 
bon consumption in dpb4A increased transcriptionally 
(9 of 12 in total), including genes for fermentation 
{IFD6), glycogen catabolism, and the xylose catabolic 
gene XYL2, The genes of these three metabolic processes 
also were upregulated in RBFl and HFLl mutants. 
Therefore, we assume that the growth defects of RBFl 
and HFLl mutants were also contributed by their re- 
duced abflity to use non-glucose carbon sources includ- 
ing lipids mentioned above. However, gene transcription 
of glycolysis and fermentation was upregulated in each 
mutant. 

Amino acid metabolism is regulated by each TR 

Regarding genes of amino acid biosynthesis, more genes 
were downregulated than upregulated for each of the 
TRKO mutants (Table 4). However, for the hfllA and 
dpb4A, down regulation of methionine synthesis genes 
were particularly common. Interestingly, transcription of 
the aromatic amino acid catabolic genes AR09 and 
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Table 4 The transcription profiles of alternative carbon utilization and phenotype-related genes among TRKOs 



Biological processes 



rbflA 



hfllA 



dpb4A 



Lipid metabolism 



Non-glucose and glucose utilization 



Amino acid metabolism 



Morphogenesis 



Transporters 



(n = 62 f 

Dw-Peroxins (4/4)^ 
Dw-lipid catabolism(29/31) 
& glyoxylate cycle (2/2) 
Dw-PL biosynthesis (10/12) 
Up-PL catabolism (3/4) 
Dw-SL biosynthesis (3/4) 
Dw-ERG biosynthesis (2/4) 
(n = 31) 

Dw-carbon utilization (26) 

GALh GALIO 
Up-fermentation 

glycolysis 

glycogen 

glucose utilization 

xylose 
(n = 31) 

Dw-aa biosynthesis (8) 
Up-aa biosynthesis(3) 
Dw-aa catabolism (5) 

Up-aa catabolism(9) 
ARO9ARO10 

Up-sulfur/nitrogen assimilation (6) 

(n = 27) 

Up-hyphal formation (13) 
ECEJ HWPiPEFh HGC1,FGR43 
RBRl, IHD2fGR6-lA,10 
(n = 101) 

Dw: sugar, amino acid, MSF 

sterol/PL, nucleosides, 

choline, nicotinamide, 

ion (K+, NHJ, Ca+^ P", CI") 
Up: urea, allantoate 

spermidine/polyamine 

cation (H^ Cu^^ Fe^^) 



(n = 52) 

Dw-Peroxins (5/5) 
Dw-lipid catabolism(14/17) 
& glyoxylate cycle(2/2) 
Dw-PL biosynthesis (15/17) 
Up-PL catabolism (2/2) 
Dw-SL biosynthesis (2/2) 
Dw-ERG biosynthesis (3/4) 
(n = 32) 

Dw-carbon utilization (23) 

GALl 10,102,7 
Up-fermentation 

glycolysis 

glycogen 

glucose utilization 

trehalose 
(n = 28) 

Dw-aa biosynthesis (9) 

MET2,3,10,15 
Dw-aa catabolism (5) 

AR03, AR08 
Up-aa catabolism(8) 

ARO9,ARO10 
Dw-sulfur/nitrogen assimilation (6) 
(n = 33) 

Up-hyphal formation (12) 
ECEl, HWPl, FGR18 ,HGC1 
FGR43, RBR1JHD2 
(n = 80) 

Dw: sugar, amino acid,MSF 
sterol/PL, nicotinamide, 
CDRs efflux pump, urea 
ion (S", NHJ, Zn+^ p-) 

Up:spermidine/polyamine 
cation (H^ Ca+^Cu^^, Fe+^) 



(n = 19) 

Up-lipid catabolism (6/9) 
& glyoxylate cycle(2/2) 
Dw-PL biosynthesis (4/4) 
Up-PL catabolism (3/3) 

Up-ERG biosynthesis (2/2) 
(n = 12) 

Up-carbon utilization (9) 

Up-fermentation 

glycolysis 
glycogen 
glucose utilization 
xylose 
(n = 19) 

Dw-aa biosynthesis (8) 
MET2,3,6,1 0,1 3,14 
Dw-aa catabolism (5) 

Up-aa catabolism (6) 



(n = 17) 

Up-hyphal formation (8) 
FGR6- 1,3,4, 10, RBRl, IHD2 

(n = 37) 

Dw: lactate, polyamine 



Up: glucose, acetate, MSF 
fatty acid, aa, 
ions (H^ Cu^^ Fe^' , S") 



a: Total number of genes in this group; 

b: x/y indicates "x" number of genes are down (Dw) or up (Up) regulated among total of "Y" number of genes in this metabolic process. 



AROlO were up-regulated only in rbflA and hfllA 
(Table 4). Both gene products are aromatic transami- 
nases [31]. Their functions are associated with providing 
an alternative, energy efficient means for NADH regen- 
eration, nitrogen assimilation, and pseudohyphal growth 
[31]. As stated above, down regulation of the MET genes 



was observed in hfllA and dpb4A, Methionine, as a 
constituent of proteins, is also critical to biochemical 
pathways, including the "methyl cycle" which generates the 
key metabolite S-adnosylmethioinine (AdoMet) [32]. As 
the main donor of methyl groups in methylation reactions, 
AdoMet plays a vital role in de novo phosphatidylcholine 
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(PC) synthesis that requires three AdoMet-dependent 
methylation steps [33]. 

Morphogenesis and cell wall responses are regulated 
by each TF 

The repressive activity of RBFl on filamentous growth 
in C albicans was first noted by Aoki et al [22]. In 
Table 4, we list the most common genes that are related 
to filamentous growth and their expression level in each 
mutant. We show that the production of hyphae was 
associated with the upregulation of genes, such as 
RBRl, HWPl and ECEl in rbfl^ and hfll^ mutants, but 
much less so in dpb4A, Transcriptional changes were 
not noted in the transcription factors CPHl and EFGl, 
These partial transcriptional profiles mostly correspond 
to the hyphal phenotypes of the rbflA and hfllA men- 
tioned above. 

Microarray data support a general increase of genes 
encoding cell wall |3-glucan biosynthesis among three 
mutants, such as EXG2, PHRl, PHR2, GSCl and KREL 
Up or down regulation of genes associated with the 
regulation of mannosylation are noted in the hfllA and 
rbflA (Additional file 1: Table SI and Additional file 2: 
Table S3). In addition to the cell wall glucan biosynthesis 
genes, those of the cell wall integrity and MAPK pathways 
were up-regulated, including the CHKl histidine kinase 
and the CEKl MAP kinase. Both genes are known to regu- 
late cell wall polysaccharide synthesis [34,35]. 

Regulation of metabolic flux transporters 

The regulatory roles of the three TRs on transporter ac- 
tivity have been noted (Tables 3 and 4). The major 
changes in both rbflA and hfllA mutants were down- 
regulation of transporters for sugar, lipid, amino acids, 
as well as the MFS transporter family (major facilitating 
superfamily). Quantitatively, 101 transporters were down- 
regulated in rbflA, 80 in hfllA, and 37 in dpb4A, of which 
the mitochondrial transporters and inter organelle trans- 
porters are not included. Certainly, the circuits for nutrient 
import from extracellular environment or intracellular 
translocation between compartments are regulated by all 
TRs but less so by DPB4, In dpb4A, gene expression for 
MFS, sugar, lipid and amino acid importers are increased. 

The measurement of intracellular accumulation of 
R6G is a useful method to reflect the activity of the CDR 
drug efflux pumps. The extracellular release of R6G in 
C. albicans was inversely correlated with the level of this 
group of efflux exporters [36]. Similar to goal A, the 
CDR genes {CDR2, CDR4 and CDRll) are down regu- 
lated in hfllA, which may explain its poor extracellular 
efflux rate of R6G shown in Figure 7 and hypersuscepti- 
bility to fluconazole (Table 2). However, these CDR 
genes were not changed in rbflA and dpb4A although 
they displayed a similar rate of R6G efflux as hfllA, 



Rbflp, Hfllp, or Dpb4 may regulate efflux by a different 
mechanism. Because R6G has a permanent positive 
charge, its cellular accumulation relies on a plasma 
membrane potential that is localized mainly in the mito- 
chondria [37]. 

The spermidine transporter was only upregulated in 
rbflA and hfllA, These data may illustrate that these 
mutants have a high demand for sustaining intracellular 
pH and membrane potential since the spermidine trans- 
porter synchronizes Ca^^, Na^, K"^ -ATPase in plant 
cells [38]. However, transporters of metal cations were 
upregulated in each of the TRKO mutants. The signifi- 
cance of uptake of Fe^^ and Cu^"^ uptake is related to 



8000 -, 




B 



E 

3 



70 n 



30' 40' 50' 60' 
Time (minutes) 

UmtDNA/nDNA 
DmumtDNA /wtnDNA 




WT dpb4A hfHA rbf1A 

Figure 7 Membrane transport of R6G is reduced in each TF 
mutant and relative mtDNA copy number is less in c/p64A 
mutant (7B). (A) Parental (SN250) and each mutant were assayed 
for transport of R6G. Relative RFU, relative fluorescent units over a 
90 min time interval were determined. Cells were starved in buffer 
for 20 min then glucose was added to each culture. Compared to 
parental cells, all mutants had little transporter activity. (B) The ratio 
of mtDNA copy number to that of nuclear DNA (nDNA) is calculated 
by ACt with matched pairs of mtDNA/nDNA primers. The relative 
copy number of mtDNA is averaged from three biological replicates 
(mtDNA/nDNA). Compared to the parental strain, dpb4A has less 
mtDNA copies compared to its own nDNA or parental nDNA. 
However mutants (rbflA and /if//A) have a similar mtDNA copy 
number as the parental strain. (mumtDNA: mutant mtDNA; wtnDNA: 
wild type strain nDNA). 
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mitochondrial respiration since electron transfer among 
ETC complexes is carried out by reduced metal ions [39]. 
The high demand for metal ion uptake in the TRKO mu- 
tants again suggests their defective mitochondria. 

Dpb4p is required for mitochondrial genome maintenance 
in C. albicans 

Giving the fact that a few complex I genes are down- 
regulated transcriptionally in dpb4A described above 
(Figure 4), we performed real-time PCRs to determine if 
mtDNA maintenance is affected in this mutant with four 
sets of primers: two sets of mtDNA encoded genes 
NADl (complex I subunit) and COXl (complex IV sub- 
unit), and two sets of nDNA genes (18S rRNA and 
SODl). The average number of copies of mtDNA per 
nDNA for DPB4 TRKO strain is less than half the levels 
of WT and other two mutants tested in Figure 7B. Since 
nDNA replication is also extensively affected in dpb4A 
mutant microarray data, we also normalized the mtDNA 



copy number by comparing the mtDNA Ct of dpb4A 
with nDNA Ct of WT cells. Again we see the reduction 
of mtDNA replication rate in this mutant. 

The TRs regulate other TRs 

In eukaryotes, -3-5% of their protein repertoire is tran- 
scription factors [40]. In our previous studies of GOAl, 
we found -100 transcription factors that were down reg- 
ulated in the gene-deleted mutant, including reduction 
of a large group of Zn2-Cys6 cluster TRs. Presumably, 
the TRs regulate metabolic pathways. However, of 77 
such genes, most are poorly characterized, yet they are 
fungal-specific [41]. We compared the regulation of 
other TRs by RBFl, HFLl and DPB4, The TRs regulated 
by RBFl and HFLl are closely related. Fifteen of the TR 
genes were either up or down regulated (Figure 8 and 
see Discussion). In this group, some genes shared be- 
tween RBFl and HFLl mutants were also changed in 
the GOAl mutant, such as ZCFl ZCFS, ZCF16, ZCF21, 




ZCF13 UME7 
ZCF15 ZCF23 
ZCF28 STP3 
ECM22 URA3 
Orf19.1757 



rbf1A 



Down- regulated pathways: 

Mitochondrial respiration 
Lipid oxidation 
Pliospliolipid biosyntliesis 
Otiier carbon metabolism 
DNA replication 
Oxidative stress response 
Hyphal morphological switch 
Cell wall response 
Adhesion 
Up-regulated pathways: 
rRNA processing 



MIG1 SUT1 
STD1 
orf19.173 
HCM1 ZCF14 
SFL2 



ZCF3 ZCF1 
ZCF16 ZCF21 
FCR1 TRY5 
Orf19.5953 
TRY4 

ZCF3 SEF2 
SFU1 CRZ1 
ETH1 TEC1 
RFX2 




WOR3 
BRG1 
CZF1 
UME6^ 
EFH1 



hfHA 



Down-regulated pathways: 

Mitochondrial respiration 
Lipid oxidation 
Phospholipid biosynthesis 
Other carbon metabolism 
Up-regulated pathways: 
DNA replication 
rRNA processing 
Oxidative stress response 
Hyphal morphological switch 
Cell wall response 
Adhesion 



Down-regulated pathways: 

Mitochondrial respiration 

Lipid oxidation 

Phospholipid biosynthesis 

Other carbon metabolism 
Up-regulated pathways: 
DNA replication 

rRNA processing 
Oxidative stress response 
Hyphal morphological switch 
Cell wall response 
Adhesion 



^pb4A 



Down-regulated pathways: 

Mitochondrial respiration 
Phospholipid biosynthesis 
DNA/mtDNA replication 
rRNA/mt-rRNA processing 
Up-regulated pathways: i 
Lipid oxidation \ 
Other carbon metabolism 
Oxidative stress response 
Hyphal morphological switch 
Cell wall response 
Adhesion ^ 



Figure 8 TRKOs and their influence on transcription of other TR genes. Each TRKO mutant and goal A is shown with arrows that connect 
common TR genes (rectangular boxes) affected in each null mutant. The TR genes are indicated in red (upregulated) or green (down regulated). 
Thus for dpb4/\, 5 upregulated genes are shown (connected by a right-facing black arrow) that are common to both rbflA and hfllA. Common TR 
genes of rbflA and hfllA are similarly inked by black arrows. The number of TR genes common to both rbflA and hfllA is much greater than 
those common to dpb4/\. Also shown are comparisons of presumed up or down regulated pathways for each mutant including goal A (square 
boxes of various colors). 
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FCRl TRY4 and RFX2. The likely involvement of these 
7 TRs in the regulation of metabolic process needs to be 
determined, although the last three transcriptional regu- 
lators have been reported to process filamentous growth 
and azole drug resistance [42,43]. 

Only 5 common gene changes were noted in the three 
TRKOs (Figure 8). Interestingly, the predicated functions 
of UME6, W0R3, BRGl, CZFl and EFHl are all associ- 
ated with hyphal formation or the white -opaque switch 
that actually matches their hyphal phenotypes. However, 
these genes were down-regulated in goalA which corre- 
lates with the inability of that strain to form hyphae 
under hyphal inducing conditions. In conclusion, RBFl, 
HFLl and DPB4 are repressors of the five genes men- 
tioned above for morphologic switching. Overall, the 
regulatory network for DPB4 mutant is very different, 
since only 14 of TRs were changed in dpb4A compared 
to 44 for rbflA and 31 for hfllA, including the five com- 
mon TRs just mentioned above and two other TRs all of 
which relate to morphogenesis. 

Discussion 

Goalp is required (directly or indirectly) for a number of 
cell processes [16-20]. Those processes include homeostasis 
of stationary phase cells, morphogenesis, prevention of 
early apoptosis, ATP production through its regulation of 
the ETC CI, and communication between peroxisomes and 
mitochondria. The latter function ensures that both organ- 
elles coordinate pathways of energy production, including 
aerobic respiration, |3-oxidation, gluconeogenesis, and car- 
bon conservation via the glyoxylate cycle. We believe these 
are the core functions of Goalp. Deletion of GQAi, 
achieved by constructing a null mutant, causes loss of 
chronological aging, constitutive yeast growth, lethal levels 
of ROS, apoptosis, heightened susceptibility to triazole anti- 
fungal drugs, and reduced energy production from both 
mitochondria and peroxisomes. The host-pathogen out- 
comes from a lack of Goalp are avirulence, a failure to 
colonize kidney tissues, and increased killing by human 
neutrophils and macrophages [16,20]. 

The availability of a transcription regulator knockout 
(TRKO) library proved suitable for facile screening based 
upon the ability of each mutant to grow or not on YP 
medium containing glycerol as the sole carbon source. 
Six such mutants were then reduced to three whose 
gene absence caused a down regulation of GOAl indi- 
cating that each TR is a positive regulator of GOAl, 
However, we have not as yet identified gene targets of 
each TF. Thus, our supposition of direct or indirect 
regulation by each TF awaits confirmation. Our ongoing 
studies focus on chromatin immunoprecipitation (ChIP) 
to identify gene targets of each TF. 

We have not completed phenotyping profiles of the 
TRKOs except for antifungal susceptibilities and respiratory 



competency. In this regard, Rbflp and Hfllp are similar to 
Goalp in that they each are required for maintenance of 
WT levels of fluconazole susceptibility. Each of the three 
TRs is required for respiratory activities but to varying de- 
grees. For example, Rbflp is much more critical to the 
maintenance of CI activity and, consequently, low levels of 
ROS. Each of the TRs is a negative regulator of morphogen- 
esis (yeast-hyphal transition), yet Goalp would appear to be 
part of a positive regulatory circuit. Further, Goalp is crit- 
ical to host cell recognition. Loss of GOAl results in a 
downregulation of genes encoding this function (17-19). 
Oppositely, the TR mutants display upregulation of these 
genes. This difference needs to be correlated with the bind- 
ing of the TRKOs to epithelial and innate immune cells. 
The cell wall inhibitor studies may point to a role of each 
TF in compensatory changes in response to cell wall defects 
caused by inhibitors. A last look at functional annotation is 
that of transporter activity, of which each of the TRs and 
Goalp is essential (Figure 7A). 

An important advantage of the TRKO mutants is to 
identify genes that are common to specific biological 
processes. For example, SUV3 is the only common 
down- regulated gene of mitochondrial genome mainten- 
ance for all TR mutants (Table 3). SUV3 is an RNA 
helicase that is essential for mitochondrial RNA catabol- 
ism and respiratory growth. Oppositely, for the filament- 
ous morphogenesis circuitry, three genes among several 
are commonly upregulated in each TR [RBRl and IHD2 
(Table 4)]. 

To visualize the overall roles of each TR and Goalp, 
we aligned genes that were up or down regulated in each 
mutant and also shared among the TRKOs. Also, for 
each TR, we have summarized the gene responses as 
presumed functional activities (Figure 8). The three TRs 
share similar roles in metabolic pathways, stress re- 
sponses, cell wall regulation, fluconazole susceptibility, 
and morphogenesis, but transcriptional changes are 
different. As stated previously, there was a higher num- 
ber of genes committed to metabolic processes in Rbflp 
and Hfllp than in Dbp4p regarding alternative carbon 
utilization of non-favorable carbon, lipid and amino acid 
catabolism. Oppositely, metabolic regulation performed 
by Dbp4p is likely because of its role in the maintenance 
of mtDNA genome (Figure 7B). Thus, when mtDNA 
replication rate in dpb4A is compromised, the cell me- 
tabolism and respiration are also then affected. While 
Rbflp and Hfllp regulated lipid oxidation and other 
peroxisomal functions such as the glyoxylate pathway, 
Dbp4p is less critical to these processes. Likely, phospho- 
lipid homeostasis in C. albicans is also contributed by each 
TR but in different ways. In fact, a more than 100-fold re- 
duction of 1N04 suggests that Dpb4p provides a positive 
regulation of PL biosynthesis. Ino4p is the PL biosynthetic 
activator that forms a heterodimer with Ino2p to regulate 
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UASiNo-containing genes [44]. When PA levels drop in ER, 
the Opilp repressor disassociates from the nuclear/ER 
membrane, then translocates to nuclei to attenuate the 
Ino2p-Ino4p complex activity [45]. 

Oxidative stress caused by ROS is usually associated 
with aberrations of mitochondrial metabolism that can 
threaten cell survival [46]. In this regard, the levels of 
cell ROS are quite different among the TRKO mutants. 
ROS production by rbflA was much higher than in 
hfllA, while that of dpb4A was similar to WT cells. 
However, most stress response genes in rbflA and hfllA 
were similarly represented. Perhaps the higher level of 
ROS in rbflA is due to the combined loss of CI and CIV 
activities. CI accounts for ~ 50% of the oxidoreduction 
activity of the ETC, and if dysfunctional, one could ex- 
pect much more ROS as superoxide. Another major 
functional difference in the three TRs as compared to 
Goalp is that they negatively regulate genes associated 
with cell wall and cell integrity while Goalp positively 
regulates this gene cluster including the signaling trans- 
duction pathway MAPK CEKl and C/Z/O [24,25]. Again, 
direct versus indirect regulatory effects need to be sorted 
out. Also, as stated above, Goalp is required for the 
yeast to hypha phase transition along with cell surface 
GPI-anchored adhesins, while each of the three TRs re- 
press these growth requirements. Finally, we note again 
that Dbp4p positively regulates nuclear DNA/mtDNA 
replication or maintenance as well as rRNA/mt-rRNA 
processing. The assignment of genes in these groups 
quantitatively is much greater in Dbp4p than the other 
two TRs. 

Translational applications of these data are apparent. 
The TPs are fungal- specific, and therefore are of reason- 
able purpose as drug targets for antifungal drug discov- 
ery. Inhibitors likely will have severe consequences on 
cell metabolism and growth. High throughput screens 
for such inhibitors might utilize heterozygote libraries of 
TRs. Hypersensitivity of strains to compounds compared 
to a parental strain indicates a haploinsufficiency re- 
sponse (decreased fitness) suggesting a drug target. Or, 
screens of the TRKO mutants with a compound library 
may be useful. In this case, a resistance phenotype likely 
indicates the target of the compound. 

Conclusions 

Our data are the first to integrate the regulation of meta- 
bolic processes by the transcription regulators Rbflp, 
Hfllp and Dpb4p. That these TRs regulate mitochondria 
and peroxisomal metabolism also is new to fungal path- 
ogens and speaks to their regulation of organelle cross - 
talking. Equally important, each TR was previously un- 
studied except for two reports on Rbflp and its role as a 
negative regulator of morphogenesis in C. albicans. Of 
particular interest is that there is little overlap in genes 



of mitochondrial functions regulated by the 3 TRs, yet 
each is functionally important to this organelle (Table 3). 
In regard to other cell metabolic activities, such as lipid 
metabolism, non-glucose utilization, and amino acid me- 
tabolism, Rbflp and Hfllp are more similar in their gene 
profiles compared to Dpb4p (Table 4). Two of the major 
differences in gene transcription among the TRs are 
transporter genes, of which Rbflp and Hfllp clearly are 
of broader scope than Dpb4p. Second, Dpb4p is a 
regulator for mtDNA genome maintenance. These ob- 
servations provide the beginning of an integrative view 
of global cell functions by transcriptional regulators of 
fungi. 

Methods 

Strains and media 

The homozygous mutants (TRKO) of a C. albicans tran- 
scription factor rbflA, hfllA and dpb4A (orfl9.2088) 
were obtained from transcription factor (TR) library pro- 
vided by Dr. Alexander Johnsons lab [21]. All assays 
with these mutants include the parental strain C. albi- 
cans SN250 [21]. All strains were grown at 30°C in 2% 
YPD, 2% glycerol YPG medium (glycerol used as carbon 
source), minimal medium (0.67% yeast nitrogen base) 
containing 2% glucose (SD), or Spider agar plates and 
YPD at pH 9.5. 

TR deletion library screens 

The entire TRKO mutant library was initially screened 
in parallel YPD and YPG plates. For these experiments, 
TR mutants with mitochondrial defects were identified 
by their lack of growth on YP-glycerol. 3 \i\ of cells was 
added to YPD or YP-2% glycerol (YPG) plates. Plates 
were kept at 30°C incubator for 24 hours and those with 
growth defects in YPG only were identified and verified 
by a second screening with a second deletion strain of 
the same gene. Subsequent experiments were done only 
with those mutants that had defective growth in YPG. 

Morphology and generation times 

All strains were grown overnight in YPD at 30°C for all 
experiments described below. Cells were washed, diluted 
to a cell suspension of 1 x 10^/ul, and streaked on YPD, 
pH 9.5 or Spider agar media and incubated at 30°C. 
Plates were observed on day 7 and photographed. The 
morphologic switch from yeast to filamentous forms in 
10% serum at 37°C for all strains was done with the 
same growth conditions. Mutants were compared to 
SN250. Generation times for rbflA, hfllA, and dpb4A 
strains were evaluated as described [16]. All strains were 
grown in YPD media at 30°C for 20 hours and cell sus- 
pensions were adjusted to an initial cell concentration of 
OD = 0.1. Also, since the mutants were constitutively 
filamentous, 50 ml of each culture was centrifuged, and 
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cell pellets were dried, and weighed every 2 hours. 
Doubling time was determined based on the biomass for 
each strain in duplicate cultures. 

Functional mitochondrial assays 

The measurement of oxygen consumption, reactive oxi- 
dant species (ROS) levels, and mitochondrial enzymatic 
activities of each mutant {rbfl, hfll and dpb4) and 
SN250 were done as described [21]. In brief, for oxygen 
consumption experiments, each strain was inoculated 
into 100 ml of YPD (2% glucose) broth until exponential 
growth was achieved. Cells were washed twice with PBS 
and suspended into fresh YPD at a cell concentration of 
OD = 0.3. 1 ml of cells was then loaded immediately into 
the sealed respirometer chamber (Hansatech Instru- 
ments Ltd., Norfolk, England). Oxygen consumption was 
measured over 10 min and polarographically recorded 
using Oxygraph Plus software. The remaining cultures 
were centrifuged to determine cell biomass. Oxygen 
consumption is presented as nmol per min per mg cell 
dry weight. Data from three experiments were averaged. 

Intracellular ROS levels for each strain were evaluated 
by staining cells using the ROS sensitive fluorescent dye 
DCFDA (2',7'-dichlorofluorescein diacetate; Sigma). 
Since growth was filamentous, the final step in ROS 
measurement was performed using a fluorescence mi- 
croplate reader in 96-well black plates (Dynex Technolo- 
gies Inc., Chantilly, VA, USA) at Xex: 485 nm and Xem: 
530 nm. Cell suspensions were kept in the dark to 
minimize loss of fluorescent signal during the assay. Cell 
cultures for each strain were prepared in 20 ml of YPD 
using an inoculum of 5 x 10^/ml; cells were grown over- 
night at 30°C, in shake culture (200 rpm). The cell pel- 
lets from 1 ml of cultures were washed once with PBS 
and suspended to 1 ml of PBS with 50 (iM DCFDA for 
30 min at 30°C, 100 rpm. Cells were washed twice with 
PBS, and 200 \A from each strain was introduced into a 
96-well microtiter plate. Cell fluorescence in the absence 
of DCFDA was used to verify that background fluores- 
cence was similar per strain. ROS data was obtained 
from duplicate cultures, and all experiments were re- 
peated a total of 3 times. 

Enzyme activities of the mitochondrial electron transport 
chain (ETC) CI and CIV were measured spectrophotomet- 
rically following procedures described previously [17,18]. CI 
(NADH:ubiquinone oxidoreductase) and CIV (cytochrome 
c oxidase) activities are plotted from duplicated samples for 
each strain as nmol per min per (ig of mitochondrial 
protein. 

Antifungal susceptibility tests 

The susceptibility (MIC50 and MICioo) for all strains to 
fluconazole, amphotericin B (AmB) and caspofungin 
was determined using the broth microdilution method 



according to CLSI guidelines M27-A3. The range of 
drugs tested was 0.25-256 (ig/ml for fluconazole; 0.03- 
32 (ig/ml for AmB; and 0.016-16 (ig/ml for caspofungin. 
Exponentially grown cultures for each tested strain were 
diluted in RPMI-1640 to a density of 1 x 10^ CFU/ml 
and 100 \A was added to each well of 96-well plate con- 
taining 100 \A RPMI-1640 with different concentration 
of drug. All plates were incubated for 48 h at 37°C. The 
MICioo was determined as the concentration resulting 
in complete growth inhibition, and MIC50 for flucona- 
zole corresponded as an inhibition of at least 50% of 
fungal growth. 

Cell wall and ETC CI and CIV inhibitor assays 

Overnight cultures of all strains were collected and 
washed twice with PBS. The cell suspension, adjusted to 
5 X 10^ to 5 X 10^ in 10 \A PBS, was spotted onto YPD 
agar with or without inhibitors. For identifying the cell 
wall defects, 25 (ig/ml of calcofluor white (CFW) or 
Congo red (CR) was added to YPD plates. CI and CIV 
inhibitors were used at concentrations of 10 (iM rote- 
none and 10 mM KCN in YPD agar. Cultures were incu- 
bated at 30°C for 24 h and photographed. 

Rhodamine 6G (R6G) efflux 

These experiments were performed using a modified 
procedure of our earlier published data [19] using 96- 
well microtiter plates. In brief, cells were initially seeded 
into 10 ml of fresh YPD after an overnight culture. Expo- 
nentially growing cells were washed twice with PBS 
(pH 7.0, without glucose), and suspended in glucose-free 
PBS to 10^/ml for 2 hours incubation to deplete glucose. 
Rhodamine 6G was then added at a final concentration 
of 10 (iM for 20 min. Again, cells were washed and 
suspended in glucose-free PBS before introducing 2% 
glucose. At every 10 min base, 0.2 ml of cells were re- 
moved and energy-dependent efflux of R6G was mea- 
sured by monitoring the absorption at 527 nm in that 
were transferred into a black 96-well plate in triplicate, 
glucose-free controls were included in all experiment. 

Quantitative PCR analysis of Mitochondrial DNA (mtDNA) 
replication rate 

The total DNAs were isolated from SN250 strain and 
mutants using RNase to remove RNA followed by stand- 
ard phenol/chloroform extraction and ethanol precipita- 
tion. The concentration of DNAs was determined by a 
nano-spectrophotometer. The primers for analysis of 
mtDNA are NADIF (5'-TAGGTTGTGTTGCTGAAT 
GTGC) and NADIR (5'-CCAGTACCACCACCCATAA 
ATAAG), COXIF (5'-GGTGAATTACGTCTAGCTGT 
TCC) and COXIR (GCACCATCTAATAGCCCTACT 
CA). Two sets of nuclear DNA (nDNA) gene are 
18SrRNAF (5'-CGCAAGGCTGAAACTTAAAGG) and 
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ISSrRNAR (5'-AGCAGACAAATCACTCCACC), SOD 
1F(5'-GCTCCAACCACAATTTCCTG) and SODIR (5 - 
TGGATTGAAATGAGGACCAGC). The 20 [iL PGR re- 
action contains Ix iQSyBR green supermix (Bio-Rad), 
0.25 [iM of each primer, and approximately 5 ng of total 
genomic DNA for each strain. PGR conditions are 
2 min at 95°G, followed by 40 cycles of 15 s of denatur- 
ation at 95°G and 30 s of annealing at 55°G and 30s of 
extension at 60°G. The relative copy number of mito- 
chondrial DNA over the nuclear DNA was averaged 
from the threshold cycle number (GJ difference for 
each pairs of mtDNA/nDNA [47,48]. The individual ra- 
tio was determined from each sets of mtDNA/nDNA 
pairs use the calculation equation N = 2^^^ where AGt = 

tical analysis of data was conducted by the t test. 

RNA and microarray analyses 

For transcriptional profiling, RNA was obtained from the 
TRKO mutants and SN250 grown in 20-ml of 2% SD 
medium at 30°G for 5 h as previously described [17,19]. 
RNA was quantified using an RNA 6000 Nano device, and 
RNA integrity was assessed using an Agilent 2100 bioanaly- 
zer. For real time PGR measurement of GOAl and NDHSl 
transcription, overnight cultures in YPD were seeded into 
20 ml of fresh SD medium containing 2% glucose. When 
exponential growth was achieved for all strains, cells were 
collected and washed, then suspended in YPG medium for 
one hour before RNA was extracted. Approximately 800 ng 
of RNA was used to prepare cDNA. Quantitative real-time 
PGR was carried out in 20 \A of Ix iQ SyBR green Super- 
mix (Bio-Rad) containing 0.25 (iM concentration of each 
primer. The experiment was performed in triplicate using 
Bio-Rad iQ5, and the transcription level of each gene was 
normalized to C. albicans 18S rRNA levels. The 2~^^GT 
method of analysis was used to determine the fold change 
in gene transcription [17,18]. 

One-color microarray-based gene expression analysis was 
done using the Agilent low input Quick Amp Labing kit. 
The RNAs for each strain were prepared from exponential 
cells cultured in 20-ml of SG medium containing 2% glu- 
cose. cDNA was synthesized from 100 ng total RNA for 
each strain according to the manufacturers instructions. 
Hybridization was completed in a Agilent SureHyb 
hybridization chamber and scan processed with an 
Agilent SGAN G2505G Microarray Scanner System. 
The array used in this study was provided by Agilent 
Technologies (eArray, ID 037557). The total of 6101 genes 
(including 12 mitochondrial genes) was done in duplicate. 

The image files were first analyzed by Agilent Feature 
Extraction Software and cyanine 3 intensities were then 
logarithmically transformed and statistically normalized. 
The fold change for each gene was calculated by com- 
paring to wild type. In this study, we adopted the cut off 



for the parametric />-value <0.05 and fold change >2 to 
determine the significance. The entire significant genes 
list for rbflA, hfllA and dpb4A are available in the supple- 
mental material (Additional file 1: Table SI, Additional 
file 3: Table S2 and Additional file 2: Table S3). 

Availability of supporting data 

The microarray data of three TRKO strains and wild 
type SN250 have been deposited to the GEO database 
with accession number [GEO:GSE54057]. The micro- 
array data of each mutant with gene changes more than 
2- fold are included in this manuscript as additional files 
indicated below. 

Additional files 



Additional file 1: Table SI. Up and down-regulated genes list in 
r6f/A mutant. 

Additional file 2: Table S3. Up and down-regulated genes list in 
c/pM/A mutant. 

Additional file 3: Table S2. Up and down-regulated genes list in 
/if/ /A mutant. 
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